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(54) Method and apparatus for reproducing three-dimensional virtual space sound 



(57) The invention provides for a method and appa- 
ratus for obtaining acoustic characteristics of sound of 
a broad frequency range of from, for example, 0 to 20 
Khz in a relatively short time with high accuracy even 
though an inexpensive computer is used; multiple sound 
ray vectors are defined; a virtual space is defined by a 
polygonal boundary: the propagation history data of the 



vector is calculated and stored, the vector being reflect- 
ed at the boundary; and, based on the data, as for each 
of the vectors, a transient response thereof at an obser- 
vation point is added to a time-series numerical array 
and stored, the response being determined on the basis 
of the reflected vector and a velocity potential deter- 
mined at the observation point by a micro-area element 
of the vector defined on the boundary. 
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. Description 

The present invention relates to a method and apparatus for reproducing three-dimensional virtual space sound. 
In particular the present invention relates to a practicable method and apparatus for reproducing acoustic charac- 

5 teristics of sound waves which are issued from a sound source and propagated to an arbitrary po'nt in a three<limen- 
sional virtual space, the acoustic characteristic appearing at the arbitrary point. More particularly, the present invention 
relates to a method and apparatus for enabling a sound having been reproduced at a predetermined point in the virtual 
space to reappear in an existing real space. Further, the present invention relates to a sound-field synthesis method 
and apparatus for very precise reproduction of a sound field in real space. 

10 Many types of methods for simulating wave motion prepagatton such as sound waves and the like propagated in 

a three-dimensional arbitrary space have been previously researched. For example, there are many conventional wide- 
use calculuses such as the sound-ray method, virtual-image method and the like based on the orthodox geometrical 
acoustics. However, these calculuses lead to considerably large calculation errors, in that they ignore wave charac- 
teristics of sound and phase information in the calculation results. On the other hand, there are some calculuses which 

15 consider wave characteristics of a sound, such as the finite element method, the boundary element method and the 
like. However, such calculuses also have drawbacks in that it is very difficult to calculafe a transient response; and, a 
very large volume of calculation is required to obtain calculation results for audio frequency bands with a high frequen- 
cies such as 1 6 Khz, and the like. These facts make it very difficult even for current super computers to obtain calculation 
results. 

20 As for the above problems, the inventor of the present invention previously proposed an approximate calculus 

(hereinafter referred to as the approximate boundary integral method) which takes into account the wave characteristics 
by using a modification of Kirchhoff's integral equation as its basic theoretical formula, the Kirchhoff integral equation 
being known as one of integral representations of three-dimensional inhomogeneous wave equations. This approxi- 
mate boundary integral has been clearly proved to be capable of realizing a very close approximation to acoustic 

25 characteristics propagated in a three-dimensional space. On the other hand, sound waves propagated in the space 
act as waves and propagate in all directions in the space, and are absorbed by wall surfaces defining the space or 
reflected at the surfaces. The sound waves are thus reflected further, and propagate in all directions in the space. In 
order to perform an acoustic analysis of such sound waves, the approximate boundary integral may be used. However, 
even by the use of this approximate boundary integral, a very large volume of wave calculation as in infinite series is 

30 still required. 

The present invention seeks to provide a sound reproducing method and apparatus having advantages over known 
methods and apparatus. 

It is an object of the present invention to provide a processing method and apparatus for reproducing acoustic 
characteristics of sound waves propagated in a three-dimensional space, which method and apparatus: use the Kirch- 
35 hoff integral or the approximate boundary integral; are free from the above-mentioned problems and excellent in ac- 
curacy and practicable; further minimize a software's static/dynamic size; and, enable even an inexpensive computer 
to calculate with high accuracy the acoustic characteristic of sound waves having a wide frequency band, for example, 
from 0 to 20 Khz in a relatively short time, without using a high-speed and expensive electronic computer such as 
super computers and the like. 

to it is another object of the present invention to provide a method and apparatus for actually reproducing acoustic 

characteristics of sound waves in an existing real space, the sound waves being issued from a sound source, propa- 
gated in a virtual space different from the real space, and affecting a desired position of the virtual space. 

According to a first aspect of the present invention, the above objects of the present invention are accomplished 
by providing: 

45 A method for reproducing a three-dimensional virtual space sound appearing at an arbitrary point in a virtual space 

defined by a plurality of boundaries, the sound being issued from at least one sound source and propagated in the 
space, the improvement wherein: 

sound waves radiating from a wave source are represented by a plurality of sound ray vectors; of the boundaries 
intersecting with the sound ray vector, as for one, which is within a distance that the sound waves travel in a predeter- 
50 mined period of time, upon which the sound ray vector is incident, and at which the sound ray vector is reflected, a 
propagation history data of each of the sound ray vectors is stored, the propagation history data comprising the incident 
sound ray vector, the reflected sound ray vector, a total propagation distance between the wave source and one of the 
boundaries, and coordinates of the intersection at which the sound ray vector intersects with the one of boundaries; 
and acoustic characteristics of the sound appearing at the observation point are determined on the basis of both the 
5S stored propagation history data and a micro-area of the one of boundaries occupied by the sound ray vector corre- 
sponding to the propagation history data. 

According to a second aspect of the present invention, the above objects are accomplished by providing: 

A method for reproducing the three-dimensional virtual space sound, as set forth in the first aspect of the present 
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invention, wherein: 

the acoustic characteristics of the sound affecting the observation point at predetermined time intervals are added 
to a time-series numerical array corresponding to the predetermined time intervals and stored, so that a transient 
response of the sound appearing at the observation point is determined. 
5 According to a third aspect of the present invention, the above objects of the present invention are accomplished 

by providing: 

A method for reproducing the three-dimensional virtual space sound, as set forth in the second aspect of the 
present invention, wherein: 

10 a plurality of loud speakers for reproducing the virtual space sound in a real space are arranged in a system; 

in accordance with positions of the loud speakers, a plurality of sound directions in which the sound waves reach 
a listener are defined in the system; and 

the transient response of the sound is determined as to each of the plurality of sound directions thus defined. 

is According to a fourth aspect of the present invention, the above objects of the pcesent invention are accomplished 

by providing: 

A method for reproducing the three-dimensional virtual space sound, as set forth in the third aspect of the present 
invention, wherein: 

in the system for reproducing the virtual space sound in the real space, the transient response of the sound is 
20 determined as to each of combinations of (at least multiple one) sound source(s) and the sound directions thus defined. 

According to a further aspect of the present invention, the above objects of the present invention are accomplished 
by providing: 

The method for reproducing the three-dimensional virtual space sound, as set forth in the fourth aspect of the 
present invention, wherein: 

25 the transient response of the sound corresponding to the loud speakers is reproduced by the use of a sumof- 

products calculator so as to reproduce a sound field at the observation point in real space. 

Further, the present invention introduces the concept of virtual windows to realize a method and apparatus for 
reproducing sound with high accuracy, and also to realize a method and apparatus of synthesis of acoustic character- 
istics of a sound, as follows. Namely, in the method and apparatus of the present invention for reproducing three- 

30 dimensional virtual space sounds: a closed space which surrounds an observation point, or a wall surface from which 
the observation point is oppositely disposed is provided; the closed space or the wall surface is divided into a plurality 
of areas which are virtual windows; and, acoustic characteristics of the sound is determined at each of the virtual 
windows. In order to realize the virtual space sound in a real space, in the method and apparatus of the present 
invention: (a multiple of) loud speakers are disposed in positions corresponding to the virtual windows; and, the acoustic 

55 characteristics of the sound having been determined at each of the virtual windows are reproduced by each of the loud 
speaker thus disposed. In other words, the acoustic characteristic thus determined are reproduced by each of the 
corresponding loud speakers so that synthesis of a sound field at the observation point is realized, whereby the sound- 
field synthesis method and apparatus of the present invention are provided. 

The invention is described further hereinafter, by way of example only, with reference to the accompanying drawings 

40 in which: 

Fig. 1 is a flowchart illustrating an entire processing procedure embodying the present invention; 

Fig. 2 is a flowchart illustrating a processing procedure "C" of Fig. 1 for determining the propagation history and 

storing the same; 

45 Fig. 3 is a flowchart illustrating the processing procedure "D" of Fig. 1 for calculating the transient response based 

on the propagation history of each of the sound ray vectors by the use of a so-called approximate boundary integral; 
Fig. 4 is a flowchart of a processing procedure for calculating, adding and storing the transient response by the 
use of the approximate boundary integral; 

Fig. 5(A) is a perspective view of the boundary, illustrating a micro-area element and its area size defined by the 
50 sound ray vector; 

Fig. 5(B) is a plan view of the boundary, illustrating the micro-area element and its area size defined by the sound 
ray vector; 

Fig. 6 is a graph of a series of the sound ray vectors and the transient responses of the direct and the reflected 
sound both to be determined, illustrating their total potential at the observation point versus time: 
55 Fig. 7(A) is a perspective view of the virtual space comprising a potential sphere surrounding a listener, illustrating 

the concept of ideal sound-field synthesis: 

Fig. 7(B) is a partially enlarged perspective view of the potential sphere shown in Fig. 7(A); 
Fig. 7(C) is a further enlarged view of the potential sphere shown in Fig. 7(B); 
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Fig. 8 is a schematic diagram of an ideal sound-field synthesis apparatus; 

Fig. 9(A) is a perspective view of both the virtual space and the real space, illustrating the relationship therebetween, 
i.e., the concept of sound-field synthesis by the use of a so-called virtual window method; 
Fig. 9(B) is a partially enlarged view of Fig. 9(A); 

Fig. 10 is a perspective view of the real space provided with a pair of wall surfaces and a floor surface, in each of 
which surfaces the virtual window is disposed; 

Fig. 1 1 (A) is a perspective view of the virtual window using 96 loud speakers which perform as acoustic generators; 
Fig. 11(B) is a view similar to Fig. 11 (A), illustrating the sum-of-products calculator with 24 channels: 
Fig. 1 2 is a perspective view of the virtual space in the electric computer and the real space, illustrating the concept 
of connecting the virtual space to the real space through the virtual window: 

Fig. 13 a perspective view of the virtual space in the real space that is provided with 4 acoustic generators for 
realizing sound-field synthesis, illustrating the relationship between these spaces: 

Fig. 14 is a planned view of the virtual space, illustrating the concept of the sound-field synthesis using 4 acoustic 
generators; 

Fig. 1 5 a planned view of the real space, illustrating the concept for finding an integral of potential from the boundary 
surface as to each directions in which such potential affect the observation point: 
Fig. 16 Js a perspective view of a concrete example of the real space using 4 acoustic generators; 
Fig. 17 is a perspective view of the real space, illustrating the concept for finding an integral of horizontal compo- 
nents of potential from a micro-area on the boundary; 

Fig. 18 is a plan view of the virtual space, illustrating the concept of a so-called approximate sound-field synthesis 
method for enlarging a listening area using 4 acoustic generators; 

Fig. 1 9 is an enlarged plan view of the listening area shown in Fig. 1 8, illustrating the concept for finding an integral 
of the potential from the boundary surface as to each of directions in which such potential affects the observation 
point; 

Fig. 20 is a planned view of the virtual space, illustrating the concept of the sound-field synthesis method for further 
enlarging the listening area using 8 acoustic generators; 

Fig. 21 is a block diagram illustrating the flow of the processing procedure for realizing the synthesis of a virtual 
space in real space; 

Fig. 22 is a view of the concept of calculation illustrated in side and planned views, by the use of which calculation 
the transient response in each of directions is calculated in the case that 4 loud speakers are used; 
Fig. 23 is a schematic view of a virtual concert hall, illustrating the relationship between a singer and the remaining 
sound sources both illustrated in side and planned views; 

Fig. 24 is a schematic block diagram of a sound-field synthesis apparatus for enjoying so-called virtual reality 
Karaoke; 

Fig. 25(A) is a graph of the response obtained by the use of conventional orthodox calculation method called the 
virtual-image method; and 

Fig. 25(B) is a graph of the impulse response obtained in the virtual space by the use of the calculating procedure 
of the present invention. 

Fig. 26 is an arbitrary space, for illustrating the concept of the Kirchhoff integral equation and called the approxi- 
mate-boundary integral equation. 

First of all, in order to facilitate the description of the present invention, the Kirchhoff integral and a so-called 
approximate-boundary integral will be described. Then, the present invention will be clarified. 

The Kirchhoff integral equation (i.e., equation 2) is derived from three-dimensional wave equation (i.e., equation 
1 ). A modified integral equation (i.e. , equation 3) of the equation 2 is called the approximate-boundary integral equation. 

With reference to Fig. 26, the three-dimensional wave equation (i.e., equation 1) is represented as follows: 

vV4>^ } =-F(x, y,z,t) 

cr dt 

where: x, y and z are three-dimensional coordinates; C is the sound velocity: t is the time: F is a wave source; and 
<t> is the velocity potential 

where: S2 is the surface defining an arbitrary space; Q is an arbitrary point on the S2: P is the observation point 
for determining the transient response; V is an area except an area V1 in the vicinity of the observation point P; F 
is the wave source; dS 2 is a micro-area vector on the S2. 

Equation 2 is an integral of three-dimensional wave equation (i.e., equation 1) to represent the velocity potential 
at the observation point P, and is represented as follows: 
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p 4n^r c rf 4nJj, r 0 Crf c « 0 r cr at 2 

5 where: (|>p is the velocity potential at the observation point P; [] td is the delay time t-r/e; V is an arbitrary space sur- 
rounding the observation point P; S2 is the surface of the V; Q is an arbitrary point on the S2; and, r is the distance 
between the Q and the observation point P 

The modification of the integral of Kirchhoff's equation (i.e., the approximate-boundary integral equation 3) is rep- 
resented as follows: 

TO 

♦ p-^fft--! / (— l - *-f (t „I°j (c0(S (7, a) -cos (r^a)| 

d C 4k * Crr Q at C 



f I * ' — ! ( ( r t a) - — cos ( r 0 n) j 

20 rr Q C r r 0 - 2 

where: <(>p is the velocity potential at the observation point P; [] td is the delay time t-r/e; V is an arbitrary space sur- 
rounding the observation point P; S2 is the surface of the V; Q is an arbitrary point on the S2; r is the distance between 
the Q and the observation point P; and, r 0 is a total distance between the wave source and the Q. 

25 The function f(t) represents a transient signal of a sound produced at the wave source. 

Now, a processing procedure according to a first embodiment of the present invention will be described. In the 
embodiment, although the present invention is described by the use of the approximate-boundary integral equation, it 
is also possible to describe the present invention by the use of any other suitable equation, provided that such other 
equation may determine the velocity potential of a sound as in the Kirchhoff integral equation, three-dimensional wave 

30 equation and the like. 

The flowcharts shown in Figs. 1 to 4 facilitate the description of this embodiment of the present invention. 
In these flowcharts, steps of the processing procedures are denoted by sequential reference numerals starting 
from the number of T. Fig. t shows the entire processing procedure ranging from the start to the end of the procedure, 
and is divided into two particular parts: "processing procedure A" and "processing procedure B". 

35 A step 7, i.e., processing procedure "C" shown in Fig. 1 corresponds to the processing procedure shown in Fig. 

2. On the other hand, processing procedure °D" shown in Fig. 1 corresponds to the processing procedure shown in 
Fig.3. Further, processing procedure °E" shown in Fig. 3 corresponds to the processing procedure shown in Fig. 4. 

As shown in Fig. 1, in a first step 1 "initialize", performed as prerequisite for executing the processing procedure 
comprises "setting of calculation conditions", "setting of wave-source conditions", "setting of a boundary", and "setting 

40 of wave-source radiating sound ray vectors". 

Set in the "setting of calculation conditions" are: three-dimensional coordinates of the wave source; the number 
of observation points: three-dimensional coordinates of each of the observation points; the temperature and moisture 
of the air in the space; an analytic frequency (i.e., the highest frequency in the transient response to be calculated); 
the duration time, T of the transient response to be calculated; and, the like. Such setting is performed by means of 

45 an input unit through which a user inputs necessary data, or by means of ah external memory unit from which the 
necessary data is retrieved. 

In the "setting of wave-source conditions", an initial value of the sound issued from the sound source is represented 
to be a transient signal by the use of a delta approximate function and its derivative, which resembles an impulse in 
shape. Such initial value may be modified if necessary in application. 

50 Further, propagation velocity of the sound is calculated based on the temperature and moisture of the air set in 

. step 1. Based on the analytic frequency, a discrete-time interval, which corresponds to a frequency equal to, or over 
twice as high as, the analytic frequency, is set according to Shanon's sampling theorem. A size of memory area for 
storing therein the transient response to be calculated, which is used in the "processing procedure B", is determined 
on the basis of both the duration time of the transient response and the discrete-time interval. 

55 in the "setting of a boundary", information as to the boundary defining a sound field is set. This represents a 

numerical space for simulating wave-motion propagation of the sound therein, and constructed of a plurality of polygons 
in plane. In this "setting of a boundary", each of the polygons is called the boundary. Set in the "setting of a boundary" 
are: the number of the boundaries: a normal, i.e. : perpendicular line of each of the boundaries: coordinates of each 
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vertex in each of the boundaries; reflection and absorption at the boundaries; and, the like. In the following description 
of an embodiment of the present invention, in order to facilitate the description, it is assumed that each of the boundaries 
only performs total reflection and total absorption of the sound over the entire frequency band. 

In the "setting of wave-source radiating sound ray vectors", the number of sound ray vectors, N, radiating from the 
wave source are calculated. In this setting, the wave motion of the sound propagated from the wave source is simulated 
to calculate both a position and a time at which the wave motion is reflected at the boundary, so that acoustic charac- 
teristics (i.e., velocity potential ^p") of the sound at the observation point are calculated. Such simulation of the wave 
motion propagation of the sound is performed by the use of numerous vectors, each of which has the same solid angle 
and radiates from the wave source in the space. Such vectors are defined as sound ray vectors in the present invention. 
The number N of the radiating sound ray vectors is set so as to let the distance between adjacent sound ray vectors 
be equal to or smaller than 1/2, preferably 1/4 of a wave length °X U of the analytic frequency. The length between 
adjacent sound ray vectors depends on the capacity of processing units and the degree of approximation in the repro- 
duction of the sound, and, therefore may be more than 1/2 of the "A." in some cases. Vectors representing the wave 
surface of the sound propagated in the space, are generally called the sound ray vectors. Of these sound ray vectors, 
ones radiating from the wave source are especially called the radiating sound ray vectors. 

In step 3, the wave motion propagated from the wave source is simulated with the use of the sound ray vectors. 
Namely first of all, a direction vector Dn of the n'th radiating sound ray vector is calculated so as to have directions of 
N pieces of the radiating sound ray vectors form the same solid angles therebetween with respect to the wave source. 

Selected in step 4 following step 3 is the boundary B having an intersection at which it intersects the sound ray 
vectors propagated from the sound source. In step 5 following step 4, when there is no boundary intersecting with the 
sound ray vectors, calculation is performed as to a subsequent one of the radiating sound ray vectors. 

On the other hand, when there is some boundary having the intersection, the sound ray vector is reflected at the 
boundary B. In this case, it is necessary to judge whether the boundary showed its front or its back to the sound ray 
vector. 

Such judgment may be accomplished by determining whether an angle formed between a normal, i.e., perpen- 
dicular line of the boundary and the sound ray vector is within an angle range of from 0 to 180 degrees. In other words, 
in case that the normal vector of the boundary is so defined as to extend outward in a virtual space, when the angle 
formed between the normal boundary and the radiating sound ray vector is within the angle range of from 0 to 180 
degrees (i.e., when an inner product thereof is positive), the boundary is defined to be the front. Otherwise (i.e., when 
the angle is within an angle range of from 180 to 360 degrees), the boundary is defined to be the back. The reason 
why the angle is defined in a manner described above is that it is necessary for the computer to judge in operation: 
whether the sound ray vector is issued from the inside of the boundary toward the front side thereof (i.e., whether or 
not the sound ray vector is reflected at the boundary) ; or, whether the sound ray vector is issued from the outside of 
the boundary toward the back side thereof (i.e., whether or not the sound ray vector merely passes through the boundary 
without being reflected thereby). 

According to the above definitions, in step 6 following step 5, when the boundary B js to the back of the sound ray 
vector, it is necessary to select any other boundary having the above-described intersection. On the other hand, when 
the boundary B is, in front of to the sound ray vector, step 6 is followed by step 7, i.e., processing procedure "C". 

In step 7, i.e., processing procedure "C", when the boundary B is set so as to completely absorb the wave motion 
of the sound in step 11, the processing procedure goes to "End", as shown in Fig. 2. Otherwise, step 11 is followed by 
a subsequent step 12 in which the total propagation distance d between the wave source and the intersection Q of the 
boundary B is calculated. Step 1 2 is followed by a subsequent step 1 3. In step 1 3, when the thus calculated propagation 
distance d is one permitting the elapsed time of the wave motion for reaching the point Q to exceed the duration time 
of transient response, step 13 is followed by "End 9 , as shown in Fig. 2. Although these steps 11, 12 and 13 may be 
interchangeable in processing order, the processing order shown in Fig. 2 is the best in efficiency to minimize the 
processing time. 

Otherwise, i.e., in step 13, when the propagation distance d does not meet the above requirement, step 13 is 
followed by a subsequent step 14 in which an incident angle "a" of the sound ray vector incident upon the boundary 
B is calculated to determine the incident angle formed between the incident sound ray vector and the normal (i.e., 
perpendicular line) of the boundary plane in the above-described equation 3. The incident sound ray vector is defined 
as one which is the n'th radiating sound ray vector reaching the boundary B. On the other hand, the incident sound 
ray vector incident upon the boundaries V times is defined as "En, rV Consequently, the "En, r" is the radiating sound 
ray vector which is issued from the wave source to reach the r'th boundary B. 

Then, step 14 is followed by a subsequent step 15. Calculated in the step 15 is the direction vector of the reflected 
one of this incident sound ray vector, the reflected one being reflected at the intersection Q. The reflected sound ray 
vector is defined to be one which is the n'th radiating sound ray vector Dn having been reflected at the boundary B 
and is found there. Of the reflected sound ray vectors, one having been reflected at the boundaries V times is defined 
as "Fn, r". Consequently, the Tn, 1 " is the n'th sound ray vector having been reflected at the boundary for the first time. 
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Step 15 is followed by a subsequent step 16. Stored in the main memory or an external memory unit in step 16 
are: arrangement No. B of the boundary; total propagation distance "d" between the wave source and the intersection 
Q; three-dimensional coordinates of the intersection Q; the direction vector "En, r" of the incident sound ray vector 
and, the direction vector "Fn, r M of the reflected sound ray vector. A series of data calculated each time the radiating 

5 sound ray vector Dn is reflected at the boundary is defined as the propagation history of the radiating sound ray vector 
Dn. The reason why the propagation history is stored in the memory is because it is necessary to calculate the velocity 
potential of the sound ray vector, which potential affects the observation point P each time the sound ray vector is 
reflected at the arbitrary point of the boundary surface S2 (shown in the above equation 3). The total propagation 
distance "d" above described corresponds to the M r 0 D of equation 3. 

Step 16 is followed by a subsequent step 17. Selected in step 17 is another or subsequent one of the boundaries 
B having a new intersection Q at which such another one intersects with the reflected sound ray vector having the 
direction "Fn, r fl . Step 17 is followed by a subsequent step 18. When there is no boundary B having the intersection 
Q, the processing procedure goes to "End", as shown in Fig. 2. When there is any other boundary having such inter- 
section Q, the step 1 8 is followed by a subsequent step 1 9. In step 1 9, when the incident sound ray vector having the 

75 intersection Q is incident on the back side of the boundary B, it is necessary to find another boundary B having the 
intersection Q. When the incident sound ray vector having the intersection Q is incident on the front side of the thus 
found other boundary B, step 19 is followed by step 11 to repeat the above-described processing procedure. 

After completion of the processing procedure "C" with respect to the n'th radiating sound ray vector, as is clear 
from Fig. 1 , the same processing procedure is repeated for each of the subsequent radiating sound ray vectors. Namely, 

20 after completion of the processing procedures as to all the radiating sound ray vectors, the processing procedure °C" 
is followed by a subsequent processing procedure "B". 

In summary, the processing procedure "A" comprising the series of steps 1 to 9 is performed to calculate and store 
the propagation histories of all the sound ray vectors. 

As is clear from the above description, an apparatus for performing the processing procedure "A" comprises: a 

25 means for storing the initialization data; a means for calculating or determining the propagation record data; and, a 
means for storing the propagation record data. 

On the other hand, as is clear from Fig. 3, calculation performed in the subsequent processing procedure tt D" is 
based on the propagation histories of n pieces of the radiating sound ray vectors, as follows: 

Namely, the processing procedure "D" starts with step 21 . In step 21 , coordinates of the first observation point P 

30 for calculating the transient response are retrieved from memory and the like. First, data as to the propagation history 
of the radiating sound ray vector of D n = 0" is retrieved in step 22 following step 21 . As for a first one of the boundaries 
B in the thus retrieved propagation history, when this one has its back side faced to the observation point R subsequent 
propagation histories retrieved in step 32 to continue the processing procedure. On the other hand, when such boundary 
B has its front side faced to the observation point R the processing procedure goes to step 25 in which a direction 

35 vector R extending from the observation point P toward the intersection Q recorded in the propagation history is cal- 
culated. Step 25 is followed by a subsequent step 26 in which the distance RD between the observation point P and 
the intersection Q is calculated. Step 26 is followed by a subsequent step 27 in which when a straight line connecting 
the observation point P with the intersection Q intersects with the remaining one of the boundaries, it is judged that 
the velocity potential from the intersection Q does not affect the observation point P As a result, the processing pro- 

40 cedure goes to step 32. On the other hand, in step 27, when the straight line connecting the observation point P with 
the intersection Q does not intersect with any other boundaries, step 27 is followed by a subsequent step 28, and 
therein it is judged whether the wave motion reaches the observation point P within a period of time T. When it is judged 
that the wave motion reaches the observation point P within the period of time T, then the velocity potential on the 
boundary B affects the observation point R and therefore step 28 is followed by the processing procedure "E\ Inci- 

45 dentally, the above judgment is formed depending on whether the wave motion travels a predetermined distance within 
the duration time T of the initialized transient response, the predetermined distance being the sum of the total propa- 
gation distance d (i.e., that between the wave source and the intersection Q) and the distance RD. 

- As shown in Fig. 4, in the processing procedure "E": the velocity potential, which is formed by a so-called micro- 
area element of the boundary represented by the intersection Q and affects the observation point R is calculated by 

50 the use of the equation 3; and, based on a time at which the velocity potential affects the observation point P, the 
transient response is stored in array. When some data has been already stored in the same location of array, the 
transient response is added to the data. More particularly, the micro-area element is so defined as to be an area formed 
in the intersection Q by a solid angle used in the definition of the sound ray vector. In Fig. 5(A), a plurality of the sound 
ray vectors are issued toward the boundary, and one (shown in solid line) of such sound ray vectors forms the micro- 

55 area element on the boundary. In Fig. 5(B), is a plan view of the micro-area element formed on the boundary. As is 
clear from these drawings, the micro-area element on the boundary varies in size and determined by an angle formed 
between the sound ray vector and the boundary plane. More particularly in data processing, the area size of the micro- 
area element on the boundary may be equal to that of a bottom area of a cone defined by both the total distance d 
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from the wave source to the intersection Q and the solid angle described above. In this case, although approximation 
becomes poor in accuracy, it is sufficient in practice. Further, in this case, since the area size of the micro-area element 
can be determined on the basis of only both the distance d and the solid angle, it is possible to increase the operation 
speed in data processing. Incidentally, this area size of the micro-area element corresponds to the term "dS2" of the 
equation 3. 

As is clear from Fig. 4, the processing procedure starts with step 41 in which the first term of the integration term 
of equation 3 is calculated. Step 41 is followed by a subsequent step 42 in which the second term of the integration 
term of equation 3 is calculated. In calculation in the embodiment of the present invention, since the reflection and the 
absorption of a sound in the boundary are defined to be the total reflection and the total absorption only the function 
f (t) used in the calculation may be defined as to be equal to the initialized transient function of the wave source. When 
each of the reflection and the absorption are a partial one, the function f(t) corresponding to the properties of the 
boundary may be determined as to the propagation history Tn, r* each time the sound ray vector is reflected at the 
boundary. 

Step 42 is followed by a subsequent step 43 in which the area size of the micro-area on the boundary is calculated 
to determine the integral approximation. 

Step 43 is followed by a subsequent step 44 in which: a product of the calculation result obtained in step 41 and 
the area size of the micro-area element on the boundary obtained in step 43 is determined; and, the thus determined 
product and the initial value of the wave source are processed by the use of a calculus of convolution transformation. 

Step 44 is followed by a subsequent step 45 in which: a product of the calculation result obtained in step 42 and 
the area size of the micro-area element on the boundary obtained in step 43 is determined; and, the thus determined 
product and the derivative value of the wave source's initial value are processed by the use of the calculus of convolution 
transformation. 

In steps 46 to 48 following the step 45, the transient response base on the integral result of the approximate 
boundary is stored in the memory and the like. In order to reproduce the acoustic characteristics of the sound at the 
observation point R there is provided a numerical array in which the transient response to be calculated is stored: Such 
array may be provided in the initialization performed in the processing procedure "A". In the array, a location j in the 
array in which the transient response is stored is determined depending on a time Dt at which the wave motion reaches 
the observation point P. Consequently, the result of the transient response corresponding to the time Dt at the obser- 
vation point P is added and stored in the corresponding location j in the numerical array. More particularly, in the step 
46, the time Dt at which the wave motion reaches the observation point P is determined based on the sum of the total 
distance d between the wave source and the intersection Q of the sound ray vector and the distance RD. 

Step 46 is followed by a step 47 in which the location j in the numerical array corresponding to the time Dt is 
determined. Step 47 is followed by step 48 in which the data in time series obtained in steps 44, 45 is added and stored 
in the corresponding location j in the numerical array. In the above addition and storing, when some other data is 
already stored in location j, the data is added to such some other data and the resultant is stored in the location j, 
whereby the same effect as that in the approximate integration is obtained. Further, by arranging the array location in 
a manner of time series, it is possible to retrieve the sound-field characteristics at the observation point P in time series 
(i.e., in the order of the array location), which enables the sound to be reproduced, whereby the method and apparatus 
for efficiently reproducing the sound are realized. 

After completion of the processing procedure "E 8 , data processing continues for the propagation history of a sub- 
sequent radiating sound ray vector When all the data of the propagation history as to the nth radiating sound ray 
vector is calculated, the same operation is conducted as to the (n+1)1h radiating sound ray vector. Then, after com- 
pletion of calculation as to all the radiating sound ray vectors, data processing for a subsequent observation point P 
is conducted in a manner as shown in Fig. 3. 

In summary, in processing procedure "D": the velocity potential of the sound ray vector affecting the observation 
point P is calculated based on the propagation history of each of the sound ray vectors by the use of integration, the 
sound ray vector being reflected at the boundary; and, the transient response at the observation point P is determined 
and stored. 

As is clear from the above, a means for accomplishing the processing procedure "D" may be any means, provided 
that the means processes the stored propagation history by the use of approximate integration. In other words, the 
means may be. a processing unit operated by the use of software comprising a series of steps; or a computer; or some 
other unit provided with hardware for executing an appropriate processing procedure by the use of the approximate- 
boundary integration. 

After completion of the processing procedure "D", step 49 shown in Fig. 1 follows it. In step 49, in order to determine 
the acoustic characteristics of a direct sound issued from the sound source and affecting the observation point P, 
transient characteristics are added and stored in the location j of the array, the location j corresponding to a time at 
which the direct sound reaches the observation point P. 

Fig. 6 shows a graph of a series of the sound ray vectors and the transient responses of the direct and the reflected 
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sound both to be determined, illustrating their total potentials at the observation point P versus time. More particularly, 
,in the processing procedure "B", the potential affecting the observation point P is calculated each time each of the 
sound ray vectors is reflected at the boundary and stored in time series so that the transient response at the observation 
point P is determined as a whole. The processing procedure "D" and adding of the direct sound are easily understood 
s with reference to the graph shown in Fig. 6. 

After completion of the entire processing procedures and determination of the transient response of a "delta" 
approximation function given as the initial value, when the result is required in the form of impulse responses, it is 
possible to obtain such impulse response by executing an inverse sum-of-products calculation in "delta" approximation. 

Next, a method and apparatus according to a second embodiment of the present invention is described and relates 
10 to the sound-field synthesis of a virtual sound of the present invention for reproducing a three-dimensional sound field 
of a virtual space in an actual space (hereinafter referred to as the real space) will be described. Incidentally, a synthesis 
sound field realized in the real space is hereinafter referred to as the synthesis sound field. 

First, an ideal sound field synthesis system will be. described. 

In order to give a listener a so-called virtual reality or the feeling that the listener feels as if he/she were in the 
is virtual space, it is ideal that all possible sound waves reaching the listener in the virtual space are synthetically produced 
in the real space. For example, as shown in Fig. 7(A), in the case that the listener is in a virtual space, the sound waves 
radiating from the wave source form the direct sounds and the reflected sounds to reach the listener all possible di- 
rections. An ideal method for synthetically producing the sound waves (which reach the listener in the virtual space) 
in the real space is as follows. 

20 Namely, first of all, as shown in Fig. 7(A), a spherical surface surrounding the listener is given, which results in the 

fact that all the sound waves reaching the listener pass through the spherical surface without fail. Also in the actual 
world, if it is possible to realize a spherical surface through which the sound waves pass, it is possible for the listener 
to feel with high accuracy the virtual space in the real space. Such spherical surface in the real space is hereinafter 
referred to as the ideal sound-field synthesis apparatus. 

2S The ideal sound-field synthesis apparatus is provided with a infinite number of micro acoustic generators in its 

surface. The micro acoustic generators continuously spread out over the entire surface of the apparatus to radiate the 
sound waves which are composed and reach the listener. The micro acoustic generator does not reflect sound waves 
issued from the other micro acoustic generators. 

It is necessary to know when and what kind of sound waves are radiated from the micro acoustic generators at a 

30 time when the wave source radiates the sound. In this case, it is necessary for an electronic computer to calculate the 
impulse response at a time when the sound waves are issued from the wave source in a condition in which the micro 
acoustic generators on the spherical surface serve as the observation points. In this connection, since it is not possible 
to treat the infinite number of the micro acoustic generators in calculation, the number of the micro acoustic generators 
is defined to be a large number "M" for convenience of calculation. 

35 As for each of the "M" pieces of the impulse responses thus calculated, an arbitrary acoustic signal is subjected 

to a convolution transformation in a sum-of-products calculator. As a result, the sound waves are radiated from the 
"M B pieces of the micro acoustic generators. The sound waves thus radiated from the "M" pieces of the micro acoustic 
generators are composed to produce a sound field of the virtual space which is defined by the spherical surface so as 
to surround the listener. 

40 As shown in Fig. 7(B), the spherical surface surrounding the listener is divided into "M" pieces of micro areas each 

of which forms an acoustic generator Fig. 7(C) illustrates the concept of the ideal sound-filed synthesis apparatus of 
the present invention, in which each of the micro areas for producing "M" pieces of sounds is replaced with a circle 
representing a loud speaker. The sound waves radiated from each of the loud speakers represented by the circles 
show in Fig. 7(C) are composed to form a wave front reaching the listener. Based on Huygens* principle, time interval 

45 at which the individual loud speakers are arranged must be within a quarter of a wave length of the highest frequency. 
For example, in case that a sound field comprising of frequencies up to 20 KHz is composed, since the wave length 
of 20 KHz is approximately 1 .7 cm, it is necessary to define the interval of the loud speakers so as to be approximately 
0.425 cm, which is however not possible in practice. As for the number of acoustic generators, the more the number 
increases, the more the continuity of acoustic generators is improved in the space. On the other hand, from an eco- 

50 nomical point of view, it is desirable to use the least possible number of acoustic generators in the apparatus of the 
present invention. 

Shown in Fig. 8 is a sound-field synthesis system which uses the largest possible number "M" of acoustic gener- 
ators. In this case, it is possible for the listener to enjoy an ideal virtual reality space produced with high accuracy. 
However, as is clear from the drawings, it is necessary to provide an equal number of amplifiers and an equal number 
55 of channels of the sum-of-products calculator to that of the acoustic generators having been arranged so as to surround 
the listener which considerably increases the cost. 

In order to solve the the above economical and technical problems, the following measures will be adopted to 
realize a more realistic sound-field synthesis according to a third embodiment of the present invention. 
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Shown in Figs. 9(A) and 9(B) is a sound-field synthesis method using the concept of a so-called virtual window. A 
simple chamber or room called a real sound field is provided, in which the sound field of the virtual space is synthetically 
produced. In case that the virtual space of a sandy beach formed in a computer is synthetically produced in the real 
space, as shown in Fig. 9(A), there is provided a window (i.e., virtual window) in a wall surface of the real sound field, 
5 through which window the virtual space is realized in the real sound field. Such virtual window is constructed of acoustic 
generators of the ideal sound field synthesis apparatus of the prelent invention shown in Fig. 7(C), the acoustic gen- 
erators being disposed in a plane of the wall surface of the real sound field. 

In this sound-field synthesis method, the sound waves propagated from the virtual space enter the real sound field 
through the wall surface thereof. As a result, the listener in the room with such window extending over the entire wall 
10 surface of the room may feels as if he/she were at a sandy beach. As is clear from Fig. 9(B), one of features of the 
sound field synthesis method of the present invention lies in the fact that the listener may assume any position in the 
room. 

On the other hand, it is necessary to calculate the impulse response for any one of the acoustic generators disposed 
in the virtual window of the virtual space. In case that the virtual space is set as if it were a sandy beach, multiple sound 
sources are used in calculation of the impulse response. For example, as for a so-caHeo* surf sound, it is necessary to 
provide multiple of sound wave sources at the beach, and calculate the impulse response as to both such wave sources 
and ail acoustic generators. Since the number of the thus calculated impulse responses is the same as that of the 
wave sources for any one of the acoustic generators in the virtual window, these impulse responses are summed up 
and allotted to each of the acoustic generators. 
20 As described above, the virtual window, of which the acoustic generators are disposed in a plane such as the wall 

surface and the like, may be easily understood when the real sound field is set in a living space. When a picture screen 
which is acoustically transparent is provided in the vicinity of the acoustic generators disposed in the plane, it is possible 
to construct a virtual reality system with sound and picture excellent in acoustic realism. 

Further, as shown in Fig. 10, a virtual window may be provided in each of the remaining wall, floor and ceiling 
25 surfaces of the room so as to substantially realize sound-field synthesis properties of the ideal sound field synthesis 
apparatus. Also as is clear from the drawings, when the picture screen which is acoustically transparent is disposed 
in a front surface of the virtual window, it is possible to additionally provide a picture in the acoustic virtual space, which 
realizes a more improved picture system. 

Incidentally, as in the description of the ideal sound field synthesis apparatus given above, as the number of acoustic 
30 generators increases, the realism of the apparatus is improved. It is preferable to continuously dispose acoustic gen- 
erators over the entire surface. As already described above, the continuity of acoustic generators depends on the 
frequency of the sound to be reproduced. In order to realize a most preferable sound-field synthesis in theory, it is 
necessary to dispose the acoustic generators at predetermined intervals of up to a quarter of a wave length of the 
highest one of frequencies of the sound waves to be reproduced. In other words, in a system for realizing a sound- 
35 field synthesis with only low frequencies, the number of acoustic generators used therein is smaller than that of acoustic 
generators used in a system for realizing a sound-field synthesis with high frequencies. 

Shown in Fig. 11(A) is a photograph of a virtual window which comprises 96 pieces of loud speakers serving as 
sound generators. These loud speakers are disposed at intervals of approximately 19.5 cm in both horizontal and 
vertical directions. This virtual window is capable of synthetically producing sound field using frequencies of up to 
approximately 435 Hz with high accuracy. Using this virtual window, an experiment of the sound-field synthesis was 
conducted. In the virtual window, there were 96 pieces of loud speakers of which: for every four speakers, one channel 
was allotted, and, therefore a sum-of-products calculator provided with 24 channels was used together with amplifiers 
to carry out a simplified sound-field synthesis method of the present invention. Shown in Fig. 11 (B) is a photograph of 
the sum-of-products calculator provided with 24 channels. In the calculation of the sound-field synthesis of the virtual 
space, a space having a width of 14 m, a height of 10 m and a depth of 20 m was imaged so as to connect to the real 
sound field through the virtual window. Such connection between the virtual space and the real sound field through 
the virtual window is shown in Fig. 12. When the virtual window is so set as to cover all the surfaces (i.e., the number 
of which surfaces is 6) of the real sound field, the real sound field is contained within the virtual space. 

In the calculation of the sound field, since one channel was allotted every four loud speakers, 24 impulse responses 
were calculated. The base point of the calculation of the impulse responses was the center of the four speakers. The 
number of the sound sources was three. In the calculation of the impulse responses, sampling frequencies of 1 KHz, 
8 KHz and 32 KHz were set. The length of the data of the impulses was 65536 data. 

Used in the listening test were: a piano which served as the sound source disposed in an anechoic chamber; male 
and female voices in narration; a drum: a flute, and the like. 

In sound-field synthesis using the impulse response with a sampling frequency of 1 KHz, communication between 
sounds reproduced by the individual loud speakers was good. Even with the sampling frequency of 8 KHz, any con- 
siderable change in such communication was not recognized, and it was possible to clearly recognize the position of 
the piano and the like disposed in the virtual space provided with the virtual window. Especially, even when the listener 
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continuously changed his position in the real space, the position of the sound in the virtual space did not change, 
enabling the sound image to be very natural. In sound-field synthesis using impulse response with a sampling frequency 
of 32 KHz, the sound image of musical instruments such as high-hat cymbals that produce sound with high frequencies, 
identification of the position became slightly poor, and at high frequencies often became noisy This was due to the 

5 fact that synthesis of the sound waves having high frequencies fails to be continuous when each of the intervals of 
generators in the virtual window is larger than a wave length of the highest frequency of the sound to be reproduced. 
Provided as a simple measure in this case is the application of an appropriate tow-pass filter to the reproduced sound. 

(n any case, sound-field synthesis according to the method of the virtual window was capable of giving the listener 
a high quality feeling of realism. 

10 it is clear that the method and apparatus of the ideal sound-field syknthesis of the second embodiment of the 

present invention and those of the third embodiment, based on the concept of the virtual window, are all easily applicable 
to the following fields. First, they are applicable as a method and apparatus for supporting sound fields such as concert 
halls, opera houses and the like. In a concert hall, the sound may often vary in acoustic characteristics when the listener 
moves to a different place with in the correct hall. In some instances, the sound normally produced on the stage reaches 

*5 the listener or audience in a form without acoustic balance. In this case, by arranging the loud speakers in an upper 
portion of the concert hall, or in a lower portion of the concert hall, or in the wall surfaces and the like, it is possible to 
compensate the sound field in acoustic balance so as to provide the ideal acoustic characteristics according to the 
present invention, which makes it possible to substantially realize the ideal acoustic characteristics of the sound. In 
other words, by actively using such idea of supporting the sound field, it is possible for the audience to enjoy the same 

20 sound effect as that of. a famous concert hall even when they are in a hall which has poor acoustic characteristics. This 
idea is also applicable to movie theatres, rooms, public squares and the like. In these places, multiple loud speakers 
are positioned so that each loud speaker serves as each of the virtual windows, allowing sound waves passing through 
the virtual windows of the virtual space to be synthetically reproduced, which enables the audience to enjoy acoustic 
realism of the sound field of the virtual space irrespective of position. Further, the present invention is also applicable 

25 any other apparatuses, such as televisions, radios, record players, compact disk players and like apparatuses, and 
still further applicable to electronic sounds in electronic pianos and electronic musical instruments and voice-generating 
media, whereby another method and apparatus for enabling the audience to enjoy virtual reality in the sounds are 
provided according to the present invention. 

Turning now to a fourth embodiment of the present invention, it will be observed that there are used at least two 

30 acoustic generators for synthetically producing the sound field. The method and apparatus of the virtual window de- 
scribed above is a practical one for realizing the ideal sound field synthesis and may realize a sound field synthesis 
having sufficient performance. However, it is further preferable from economical point of view to reduce the number of 
acoustic generators and the number of the channels of the sumof-products calculator used for such sound field syn- 
thesis apparatus. 

35 Fig. 1 3 shows the relationship between the virtual space and the real space used in sound field synthesis using 4 

pieces of the acoustic generators. If the sound waves (which are propagated in the virtual space numerically constructed 
in the electronic computer) are synthetically produced using 4 pieces of the loud speakers disposed in the real space, 
it is possible for the listener to enjoy acoustic realism as if he is in the virtual space. 

In comparison with each of the ideal sound field synthesis apparatus and the method of sound field synthesis using 

40 the virtual window, the number of acoustic generators used in the method shown in Fig. 3 is considerably reduced, 
which makes it difficult to synthetically produce a sound field that is acoustically realistic, even when the impulse re- 
sponse of the virtual space is merely calculated at positions corresponding to those 4 acoustic generators and used 
in sound field synthesis. The reason for the difficulty described above is that it is not possible to synthetically produce 
with high accuracy the sound waves propagated between acoustic generators having been positioned apart from each 

45 other. In order to solve this problem, as shown in Figs. 1 4 and 1 5, a method for calculating the sound field for synthetically 
producing an approximate sound field by using an acoustic generator having been positioned apart from each other 
is provided. As is clear from Figs. 14 and 15, in this method, there are 4 acoustic generators. Incidentally, each of Figs. 
14 and 1 5 is a plan view of the observation point. 

Since sound field synthesis is realized using 4 acoustic generators, the potentials from the boundary surface reach 

50 the observation point (i.e., the listener's position) and are integrated in each of four directions so that 4 impulse re- 
sponses are calculated. Such integration is conducted in a manner shown in Fig. 15. In KirchhofT's integral or in the 
approximate-boundary integral, the potential from the micro area on the boundary, which reaches the observation point, 
is integrated. Namely, as shown in the drawings, the potential from the micro area affects the observation point. In 
case that the potential is propagated between the acoustic generators 1 and 2, that potential is divided into two parts 

55 in the following manner, each of which part is integrated as to each of the acoustic generators 1 and 2. As shown in 
Fig. 15, the above dividing is conducted using angles "a" and "|5 B of the observation point. As for the impulse response 
to be calculated for acoustic generator 1 , a ratio of "|5 / (« + |5)" of the potential is integrated. On the other hand, lor the 
impulse response to be calculated for acoustic generator 2. a ratio of "a / (a + |5)" of the potential is integrated. The 
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same calculation is conducted as to any other potential propagated between the remaining acoustic generators, where- 
by 4 pieces of the final impulse responses are obtained. Also, as for the direct sound, the same calculation as is in the 
above is conducted. 

In the method shown in Figs. 14 and 15, since 4 acoustic generators are disposed in a horizontal plane, it is 
essentially not possible to compose the sound waves in vertical directions. Under such circumstances, calculation is 
conducted by assuming that all the sound waves issued from a wave source (which is out of the above horizontal 
plane) to the observation point reach the observation point in horizontal directions. As a result of a listening test that 
is conducted, a sound source such as a piano and the like is in the horizontal plane, it is possible to for the listener to 
sufficiently enjoy acoustic realism without feeling any undesirable effect. Shown in Fig. 16 is photograph in which 4 
acoustic generators are disposed in a fixed manner. 

Incidentally, as -shown in Fig. 17, in case a vertical angle "0" is formed between the horizontal plane and a line 
connecting the micro area of the boundary surface with the observation point, i.e., in case that the potential with the 
vertical angle "6" reaches the observation point, a horizontal component of such potential may be calculated using a 
function of cose, and the thus calculated horizontal component may be integrated. 

In order to synthetically produce a natural position of the sound image and the reflected sound also in vertical 
directions, at least one additional acoustic generator is provided over the horizontal plane, i.e., over such 4 acoustic 
generators. Integral calculation in this case of the provision of the additional acoustic generator may be conducted in 
the same manner as that of the case in which only 4 acoustic generators are disposed in the horizontal plane. 

Now, an approximate sound field synthesis for enlarging the listening area of at least 4 acoustic generators is 
described and which relates to a fifth embodiment of the present invention. 

Clarified in the sound field synthesis using 4 acoustic generators described above are calculation of the sound 
field and the method of sound field synthesis in case that only one listener is in the sound field. However, in some 
applications, it is better from an economical point of view, for the sound field synthesis system to be modified so that 
the system is capable of admitting a plurality of the listeners to the sound field, which also enables the system to be 
applicable in various fields. Figs. 18 and 19 show a method of approximate sound field synthesis for enlarging the 
listening area using 4 acoustic generators. 

This method is different from the above-described method of sound field synthesis using 4 acoustic generators in 
the following point: namely, in this method shown in Figs. 18 and 1 9, the number of the observation point increases to 
4, at each of which point Kirchhoff's integral is conducted in the virtual space, the thus increased observation points 
corresponding in position to 4 acoustic generators. In other words, the position of each of the 4 sound generators in 
the real space is equal to that of each of four observation points in the virtual space. 

When the potential from the micro area of the boundary as shown in Fig. 19 affects the observation points, i.e., is 
propagated between the observation points 1 and 2, such potential is divided into two parts each of which is integrated 
for each of the observation points 1 , 2. As shown in Fig. 19, such dividing is conducted using angles V and "P" of a 
central point of the listening area surrounded by 4 acoustic generators. More particularly, for the impulse response 
calculated at observation point 1 , a ratio of "p / (a + p)" of the potential is integrated. On the other hand, for the impulse 
response calculated for observation point 2, a ratio of "a/ (a + p)" of the potential is integrated. The same calculation 
is conducted as to any other potential propagated between the remaining acoustic generators, whereby 4 final impulse 
responses are obtained. Also, for the direct sound, the same calculation as above is conducted. 

According to the method of calculation described above, sound field synthesis on a stage of a multipurpose hall 
was conducted, the stage being simple in shape and relatively small in size. The multipurpose hall had a width of 14 
m, a height of 15 m and a depth of 20 m, and the stage thereof was provided with a listening or listener's area of 3 m 
x 3 m. Since this multipurpose hall had a reverberation time of approximately one second, 65536 data was set for the 
impulse response. 

The wave sources disposed in the virtual space were a pair of proscenium loud speakers. 

When popular music and the like was reproduced from the proscenium loud speakers in the virtual space, it was 
found that the listener enjoyed the sound field in size and acoustic realism with much higher accuracy than that in the 
case of the above-described sound field synthesis in which 4 acoustic generators were used. 

Shown in Fig. 20 is an example of sound field synthesis which is realized using 8 acoustic generators. As shown 
in Fig. 20, eight observation points P1-P8 are arranged so as to form a circular shape surrounding a listening area. 
This shape may be modified to any other suitable shapes such as rectangular shapes, square shapes and the like. 
The 8 observation points P1 -P8 correspond to the acoustic generators disposed in the real space. Used in the example 
shown in Fig. 20 is the same method of calculation as that used in the case of the above-described sound field synthesis 
using 4 acoustic generators, so that the impulse response is calculated at each of the observation points. When the 
sound field synthesis was conducted using the same virtual space as above, it was found that the listener enjoyed a 
sound field that had much higher feeling of acoustic realism than that in the case of the above<iescribed sound field 
synthesis using 4 acoustic generators, and that the listener identified the position of the wave source more clearly. 

Next, a sixth embodiment of the present invention will be described, in which the present invention is applied to a 
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so-called "Karaoke". 

Karaoke is an entertainment in which a user sings his or her favorite song with accompaniment using a microphone 
connected to an echo machine, as if the user were a famous singer The user of Karaoke wants to feel as if he/she 
were a professional singer singing on the stage of a concert hall. Further, for musical instruments and songs, education 

5 and training may become more effective by the use of Karaoke, for the user can enjoy an ideal acoustic environment 
of being realized on the stage of a concert hall. On the other hand, today's Karaoke is a system for merely reproducing 
the user's song with accompaniment through an echo machine and loud speaker. and : therefore is not capable of 
providing an environment in which the user can feel as though he/she. were singing a song on the stage of a concert 
hall or in a music studio. This is because the acoustic difference, between a three-dimensional live sound field such 

10 as that created in an actual concert hall, and a plain sound field realized by the use of sound that is reproduced a loud 
speaker after passing through an echo machine, is considerably large. In the present invention, a three-dimensional 
sound field, which is substantially the same as that created in an actual concert hall, is synthetically produced around 
the user of Karaoke, enabling him/her to enjoy virtual reality, i.e. , feel as if he/she were on the stage of a famous concert 
hall and the like. Further, when a user plays a musical instrument, the acoustic characteristics of a famous concert hall 

15 can be re-created according to the present inventino and it is possible for the user.to acquire excellent musical edu- 
cation. 

The above is realized according to the present invention as follows: namely a process showing how the sound 
waves propagated in the three-dimensional sound field (which is to be synthetically produced) reach the user, is proc- 
essed by an electronic computer in a manner as previously described in the present invention (hereinafter referred to 
20 as sound field simulation); and, based on the thus obtained acoustic information, a virtual sound field is synthetically 
produced in the real space by the use of the sum-of-products calculator, digital-to-analog converter, loud speaker, 
amplifiers and the like. 

The above acoustic information is obtained by calculating the process showing how the sound waves propagated 
in the sound field (which is numerically constructed in the electronic computer, and, therefore hereinafter referred to 

25 as the virtual space or virtual sound field) reach the user. Consequently, the acoustic information is of the transient 
response such as the impulse, response. Further, in the transient response, it is necessary to represent the sound 
waves reaching the user in each of the directions. This calculation is executed based on the number and the positions 
of the loud speakers used when the virtual sound field is synthetically produced in the real space. 

The transient response is obtained by calculating the sound by means of the sum-of-products calculator, the cal- 

30 culator being provided with a plurality of channels corresponding to the number of the loud speakers and the number 
of the sound sources in the virtual space. 

The actual method of the present invention for reproducing the sound in a virtual space will be described with 
reference to a flowchart shown in Fig. 21 . 

First of all, various data is set before the computer simulation of the sound field is carried out. Setting contents in 

35 the above comprise: setting of the virtual space in step 51; numerical representation of the virtual space in step 52 
following step 51; and, determination of hardware to be used for sound field synthesis in step 53 following step 52. In 
step 51 for setting the virtual space such as various types of concert halls which the user wants to enjoy is designed 
so as to determine or set : the shape and size of the space; the number and positions of sound sources; required 
frequency bands; and, the listener's position. Registered in the computer in step 52 for numerically representing the 

40 virtual space thus designed are: the boundary such as the wall surface; objects disposed in the sound field; types and 
material of the boundary; and, the speed of the sound. These data enable the computer to calculate. In step 53 for 
determining the hardware, used in sound field synthesis, in order to have the user enjoy virtual sound in real space, 
the number and positions of the loud speakers are determined or set. The processing order of steps 51 to 53 may be 
inter-changed. 

45 in step 54 folbwing step 53, the simulation of the sound field is carried out so that the transient response of the 

sound propagated in the virtual space is calculated, the virtual space depending on the type of the sound field synthesis 
hardware. The process described in step 54 is substantially the same as that described in the preceding embodiment 
with the exception of the number and positions of the sound sources. Such process is described again in the following 
embodiment of the present invention for better understanding of the invention. 

50 in step 55 following step 54, the virtual space is synthetically produced in real space. The transient response is 

input to the sum-of-products calculator so that the sound field is synthetically produced by the use of sound field 
synthesis hardware, whereby the sound field of the virtual space is reproduced, which enables the user to experience 
a virtual sound field. 

Now, the apparatus for reproducing virtual sound, carried out by the method for reproducing virtual sound as 
55 described in the preceding embodiment of the present invention, will be described in a clear manner. For example, the 
method of the present invention described above is applicable to apparatus tor synthetically producing the sound field 
with reality, the apparatus being disposed around the userol Karaoke. 

In order to synthetically produce a sound field that is acoustically realistic it is important to consider the following 
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matters 1 to 4: namely, 

1. It is important to reproduce the sound that reaches the user, comprising of the user's singing voice with accompa- 
niment, that is reproduced and issued from electro-acoustic systems such as proscenium loud speakers, side-wall loud 
speakers, stage loud speakers and like systems normally installed in concert halls, multipurpose halls and the like. 

5 , For example, in general, the singing voice (for example, a popular song) of the user singing op a stage of a virtual 
space reaches a mixer after passing through a microphone and an echo machine, in an actual space. Then, the user's 
singing voice is mixed with the accompaniment of the musical band and the like in the mixer, amplified in the power 
amplifier, and then issued from the proscenium loud speakers, side-wall loud speakers, stage loud speakers and the 
like. At this time, the sound issued forward from the proscenium loud speaker radiates to the audience and a rear 
10 portion of the concert hall, and is reflected at the auditorium, rear walls, side walls and the ceiling of the hall, so that a 
part of the reflected sound returns to the user singing on the stage. 

In the remaining directions except one in which the proscenium loud speaker faces toward the audience, for ex- 
ample, in the direction towards the stage, the musical sound generally biased towards its lower-tone side depending 
on directional frequency, characteristics of the proscenium loud speaker are radiated together with the user's voice to 
reach the user on the stage in a relatively short time. In order to improve the sound in acoustic realism, it is important 
to simulate what the user hears concerning both, the sound which passes through the audience area and the sound 
which is issued from the loud speaker directly to the user. This is true for the remaining loud speakers such as the 
side-wall loud speakers installed in the side walls of the hall, stage loud speakers and monitor loud speakers installed 
on the stage and the like. 
20 in order to realize the above, the following steps (1), (2) and (3)are required: 

(1 ) Acoustic simulation of sounds with directional frequency characteristics of multiple sound sources such as loud 
speakers installed in virtual space is carried out to calculate the transient response in individual directions in which 
the sounds reach the user; 

25 (2) The sounds radiated from multiple sound sources are calculated in reflection, absorption, diffraction, transmis- 

sion and diffusion occurring at the boundary which forms the virtual space constructed in an electronic computer, 
the boundary comprising of, for example, walls and the audience of the concert hall t auditorium, floor, stage, ceiling, 
reflector boards, stage curtains, other types of curtains and the like; and 

(3) In the above acoustic simulation, the transient response is calculated with consideration for attenuation of the 
30 sound caused by atmospheric damping effect. 

2. It is also important to further improve the sound the user hears in acoustic reality by simulating the applause of the 
audience so as to synthetically produce the sound in the real space. As is in actual cases, communication between 
the singer and the audience is carried out through applause and the like. Consequently, it is very important to simulate 

3$ applause and the like for the Karaoke system. In some cases, it is also often necessary to synthetically produce some 
kind of noise such as chattering and the like occurring in the entire auditorium. 
The above is realized by the use of the following steps (1), (2) and (3): 

(1) Multiple sound sources are set in the auditorium in the virtual space to serve as the virtual audience. Sound 
40 propagation from the sound sources (i.e., audience) to the user on the stage is simulated, so that the transient 

response is calculated in the individual directions in which the sounds reach the user. In acoustic simulation, the 
transient response is calculated with consideration for the attenuation in sound caused by atmospheric damping 
effect; 

(2) When it is required to improve the acoustic simulation described above from an economical point of view, it is 
45 possible to reduce the number of the sound sources (i.e., audience) by the use of a stereo-type sound source 

constructed of a pair of sound sources installed in the auditorium of the virtual space. Namely, this is simulation 
of a stereo audio system for reproducing applause in the virtual space; and 
v (3) In order to lend higher acoustic realism to the sound field, the timing of the applause produced by the virtual 
audience can be automated. The beginning and the ending of the user's song are automatically detected to control 
s o the timing of the applause of the virtual audience. 3. As is in actual case, in some cases, it is necessary to simulate 

players of a musical band and the like in a position behind or in front of the user facing the auditorium on the stage. 
In general, the players are behind the user. However, in operas, they are in front of the user on the stage. 



55 



The above is realized by the following steps (1) and (2) : namely, 

(1) Multiple of sound sources corresponding to individual musical instruments are disposed in their positions so 
that the transient response is calculated in each individual directions in which the sounds reach the user In this 
calculation., the directional frequency characteristics of the musical instruments are also calculated: and 
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(2) From an economical point of view, a simple method is used in which a stereo-type loud speaker is simulated 
in an area of the players who are behind or in front of the stage in the virtual space, so that the sound issued from 
the stereo-type loud speaker represents the players and the musical band. 

5 4. Some music requires the users to perform a duet. In this case, each of the users often wants to sing together with 
a famous singer. Therefore, it is necessary to produce a virtual image of a famous singer, whose image is adjacent to 
the user on the virtual stage. 

The above is realized by the following step (1 ): namely, 

w (1) Acoustic simulation is carried out in a condition in which a sound source is installed in the position of the virtual 

singer (hereinafter referred to as the duet singer) adjacent to the user, so that the transient response is calculated 
in individual directions in which the sounds reach the user. At this time, a part of the duet singer's singing voice is 
also issued from another simulated sound source such as loud speakers and the like. Now, sound field synthesis, 
which is carried out using the 4 loud speakers according to the above method will be described. As for the acoustic 

is simulation, approximate calculation is carried out based on the approximate boundary integral represented by 

equation 3 which is obtained based on the assumption that the sound source is a point source. N&riables used in 
the equation 3 are the same as those shown in equation 1 . 

In case that there are, for example, 4 loud speakers (i.e., acoustic generators) which are used in synthetic pro- 
20 duction of the sound in the real space, as shown in Fig. 22, the transient response such as the impulse response of 
sound issued from the sound source to the point P is calculated in each of four directions. At this time, a sound field 
with directional frequency characteristics of the sound source is simulated. In the case 8 loud speakers are used for 
synthetically producing the sound in real space, the transient response such as the impulse response is calculated in 
each of eight directions in substantially the same manner as above. In case that the number of the loud speakers is 
25 not the same as above,- the transient response is determined in substantially the same calculation manner Further, it 
is also possible to obtain better effects by setting loud speakers in three-dimensional arrangements. 

In order to realize close simulation of what the user hears, as to the indirect and the direct sound produced by the 
electro-acoustic system such as various types of loud speakers and the like installed in the virtual space, it is necessary 
to determine how the sound waves (which are radiated from the loud speakers) reach the user on the stage by carrying 
30 out calculation with proper consideration for directional frequency characteristics of the loud speakers. Further, for the 
applause of the audience and musical instruments on the stage, it is preferable to consider their directional frequency 
characteristics. Still, further, in case that the user plays a musical instrument, the directional frequency characteristics 
of such musical instruments need to be considered. 

As for the directional frequency characteristics of loud speakers and of musical instruments, it is a great conven- 
es ience to use measured values thereof. The measured values are obtained at measuring points which are disposed in 
a three-dimensional arrangement around the sound source at angular intervals of 10 degrees. The impulse responses 
measured at the measuring points are used in calculation. 

In the approximate boundary integral according to equation 3 which is an approximation of Kirchhoff's integral, a 
way of propagation of the sound waves radiated from the sound source through the virtual space is traced. The sound 
40 waves are represented by a large number of dots forming a wave front. It is considered that each of the large number 
of such dots of the wave front carries the initial value (normally, pulses) and is propagated. At this time, when the sound 
waves (i.e., dots in the wave front) are reflected from the reflecting point of the boundary defining the virtual space, 
the velocity potential from such reflecting point affects the observation point, based on which the transient response 
of the sound waves at the observation point is calculated. At this time, with respect to the initial value, the impulse 
45 response (i.e., measurements of the directional frequency characteristics of the sound source in each of the directions) 
is calculated by the sum-of -products calculator. 

When the sound waves are reflected at the wall surface of the virtual space, a part of the sound wave is absorbed 
by the material of the wall surface. In order to calculate the acoustic influence of the wall-surface material, with respect 
to a value of a representative point of the wave front, reflection characteristics of the wall-surface's material are cal- 
50 culated by a sum-of-products calculator. It is preferable to have the reflection characteristics of the wall-surface's ma- 
terial correspond to the impulse response. This impulse response corresponds to a reflected part of the sound, the 
reflected part being measured as one impulse response when an impulse is incident on the material at every angle. 
However, normally, it is impossible to measure such impulse responses at every angle with respect to various types 
of the materials in practice. Therefore, in actual use, an approximate impulse response is calculated based on the 
5$ material's sound absorption coefficient which is determined by the use of a so-called reverberation method of the sound 
absorption coefficient. 

In the embodiment of the present invention, as shown in Fig. 23, a simple shoe-box type virtual concert hall is 
used, the hall having a width of 16 ra a height of 13 m and a depth of 25 m. Loud speakers disposed in the virtual 
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space comprised of proscenium loud speakers disposed in an upper portion of the stage and monitor loud speakers 
disposed on the stage, each of which loud speakers was provided with a right and a left channel. Each of the loud 
speakers has carefully considered directional frequency characteristics. The material used for the floor and the wall 
surfaces of the hall is hard wood. The hall is provided with multiple of acoustic absorption portions for preventinq 
s acoustic trouble such as echoing. 

The sound source of the applause that was used was actually recorded. The number of positions of the sound 
source of the applause in the auditorium of the virtual space is forty Further, for the sound field synthesis apparatus 
from an economical point of view, to reduce the number of the channels of the sum-of-products calculator, there were 
two types of applauses used. 

io As for the number of channels (i.e., the number of filters) of the apparatus carrying out convolution transform 

calculus (hereinafter referred to as the convolution-transform apparatus), the lower the number, the lower the manu- 
facturing cost. Consequently, it is preferable to reduce the number of types of sound that will be radiated in the virtual 
space. In the case that the sound field is synthetically produced using 4 loud speakers, as it is in the embodiment of 
the present invention, four channels need to be processed by the convolution-transform apparatus for a single sound 

15 source. 

The sound sources installed in the virtual space in the embodiment of the present invention are as follows: 

(1 ) the right channel on the proscenium bud speaker; 

(2) the left channel on the proscenium loud speaker; 

(3) the right channel on the stage monitor loud speaker; 

(4) the left channel on the stage monitor loud speaker; 

(5) twenty sound sources type 1 for producing the applause in the auditorium; and 

(6) twenty sound sources type 2 for producing the applause in the auditorium. 

As a resutt, the total number of the sound sources calculated for the sound simulation is forty-four. Consequently, 
the number of the sound sources to be calculated for the sound simulation is forty-four. On the other hand, the total 
number of types of sound is only four, two of which are the right and left channel of Karaoke music mixed with the 
user's singing voice, and the remaining two of which are the applause. 

As for each of the 4 transient responses calculated using a sound source comprising of the right channel of the 
proscenium bud speaker and the right channel of the stage monitor loud speaker, since an audio signal prepared by 
mixing sound of the right channel of a stereo player and the like with the user's singing voice, is calculated by the use 
of convolution transform calculus, the 4 transient responses are summed up in each of the directions in which the 
sound reaches the user. 

As for the 4 transient responses calculated using the sound source comprising of the left channel of the proscenium 
loud speaker and the left channel of the stage monitor loud speaker, the same calculation as the above is carried out. 
As a result, the total number of transient responses of Karaoke music and the user's singing voice is eight. 

On the other hand, the transient response is calculated in the forty positions of the sound sources of the applause. 
For the same reason as that described above, forty-four transient responses are obtained by using sound sources 
disposed in twenty positions, and are summed up in each of the directions in which the applause reaches the user. As 
40 a result, the total number of transient responses of the applause is eight. 

Therefore, the number of the channels of the convolution-transform apparatus is sixteen. This is because a transient 
response is required in each of four directions for a single sound source. Incidentally, even when Karaoke music is 
monaural, it is preferable to process it in the same manner as that described above. 

Shown in Fig. 24 is a schematic block diagram of the sound-field synthesis apparatus used in the embodiment of 
45 the present invention. The sum-of-products calculator show in Fig. 24 comprises of a digital-to-analog converter (D/A 
converter) and an analog-to-digital converter (A/D converter). Incidentally, the number of loud speakers used in the 
synthesis apparatus is eight. In the case that a transient response is calculated in each of eight directions for each of 
the sound sources, since the number of types of the sound sources to be calculated by the use of the sum-of-products 
calculator is four, it is necessary to provide thirty-two (i.e., 4 x 8 = 32) channels in the convolution transform apparatus. 

As described above for the embodiment of the present invention, in order to reduce the manufacturing cost by 
reducing the number of convolution-transform channels of the sum-of-products calculator, it is necessary to lower the 
number of types of the sound sources. 

The above system shown in the embodiment of the present invention was applied to Karaoke. As a result, users 
of Karaoke enjoyed acoustic realism, (which conventional Karaoke systems have never realized) as though he/she 
were on a real stage. Further, in the case that a user played a musical instrument, he/she enjoyed the same acoustic 
realism as above. 

Although the present invention has been described using an application of the invention to a Karaoke system 
above : the present invention is also applicable to any other systems described above. 
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It is possible with the present invention to analyze acoustic characteristics of sound and systems within a very 
short period of time without using a large-scale computer, which makes it possible to apply simulation of sound with 
its wave characteristics for practical use. 

Shown in Fig. 25(A) is an example of a transient response of a hall, which was obtained by using the conventional 

5 orthodox calculus called the virtual-image method. On the other hand, shown in Fig. 25(B) is an example of an impulse 
response calculated by using a program prepared according to the present invention. As is clear from these drawings, 
the orthodox calculus failed to calculate negative waves, and merely showed atmospheric exponential attenuation of 
energy. On the other hand, in the present invention, it was found that: the impulse response obtained by using the 
present invention showed boundary waves represented by negative values: phase information was calculated; and, 

10 complex attenuation was entirely calculated. 

In the program based on the processing procedure of the present invention, it was possible to accomplish all the 
necessary calculation of the program within a very short period of time with the use of a small-scale computer In 
contrast with this, the conventional calculus method required tens of thousands of hours in calculation of the impulse 
response with audio frequency band even when a high-speed computer such as super computers and the like were 

is used, and, was therefore not applicable to practical use. The present invention enables the sound to be synthetically 
produced and also to be reproduced for practical use. Since the present invention is capable of simulating propagation 
of sound with its wave characteristics, it is possible for the present invention to enable users to enjoy acoustic realism 
with higher accuracy in any arbitrary space. More specifically, the present invention provides a fundamental technique 
for the following application: namely, architecture the present invention enables users to calculate and evaluate various 

20 physical-quantity data comprising acoustic characteristics of a building such as concert halls, studios, listening rooms 
and the like before the building is constructed. 

Further, the present invention enables users to simulate desired sound fields of the interior otaircraft and vehicles, 
and also to perform various types of other simulation with reality Still further, the present invention provides other 
fundamental techniques of simulation, which enables users of the present invention to research propagation of noise 

25 produced in air ports, railways, roads, factories and like installations, and make an accurate estimate of noise influence 
upon cities, buildings, the interiors thereof and the like. 

Further, the present invention provides virtual sound reproduction apparatus for synthetically producing three- 
dimensional virtual sound field in an actual or real space, virtual sound field being substantially the same as that of an 
actual concert hall and tike installations and enabling users of the present invention to enjoy virtual reality in real space 

30 as if he or she is actually on the stage of a famous concert hall and the like. Consequently, the present invention is 
applicable to: Karaoke systems; practical apparatuses for musical instruments, songs, dances and the like; and ap- 
propriate acoustic virtual reality systems. As is clear from the above, the present invention is applicable to musical 
education and may considerably improve it in quality. 

35 

Claims 



1 . A method of reproducing a three-dimensional virtual space sound at an observation point in a virtual space defined 
by a plurality of boundaries, said sound being issued from at least one sound source and propagated in said space, 

40 characterized by the steps of representing sound waves radiating from a wave source by a plurality sound ray 

vectors; 

storing propagation history data of a sound ray vector, which data relates to one of said boundaries which is 
within a distance that said sound waves travel in a predetermined period of time and upon which said sound 
4 $ ray vector is incident and at which said sound ray vector is reflected, wherein said propagation history data 

comprises said incident sound ray vector, said reflected sound ray vector, a total propagation distance between 
said wave source and said one of said boundaries, and coordinates of said intersection at which the sound 
ray vector intersects with said one of said boundaries; and 

determining acoustic-characteristics of said sound appearing at said observation point on the basis of both 
50 said stored propagation history data and a micro-area of said one of said boundaries occupied by said sound 

ray vector corresponding to said propagation history data. 

2. A method as claimed in Claim 1, wherein said acoustic characteristics of said sound affecting said observation 
point at predetermined time intervals are added to a time-series numerical array corresponding to said predeter- 

55 mined time intervals and stored, so that a transient response of said sound appearing at said observation point is 

determined. 



3. A method as claimed in Claim 2, wherein a plurality of loud speakers for reproducing said virtual space sound in 
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a real space is arranged in a system; 

a plurality of sound directions in which said sound waves reach a listener are defined in said system and in 
accordance with the positions of said loud speakers; and 

said transient response of said sound is determined for each of said plurality of sound directions thus defined. 

A method as claimed in Claim 3, wherein in said system for reproducing said virtual space sound in said real space, 
said transient response of said sound is determined for each of the combinations of said at least one sound source 
and said sound directions thus defined. 

A method as claimed in Claim 4, wherein said transient response of said sound corresponding to said loud speaker 
is reproduced by means of a sum-of-products calculator to reproduce a sound field at said observation point in 
said real space. 

A method of reproducing a three-dimensional virtual space sound, wherein sound waves are issued from at least 
one sound source and propagated in a virtual space defined by a plurality of boundaries so that said sound is 
reproduced at an arbitrary observation point in said virtual space, characterized by the steps of disposing at least 
one virtual window opposite from said observation point and which virtual window is divided into a plurality of areas; 
and 

20 determining acoustic characteristics of said sound passing through each of said areas. 

7. A method as claimed in Claim 6, wherein a plurality of loud speakers are disposed in positions corresponding to 
said areas of virtual window to reproduce said acoustic characteristics of said sound being determined for each 
area of said virtual window, whereby sound field synthesis at said observation point is realized 

25 

8. A method as claimed in Claim 6, and including the steps of representing sound waves radiating from at least one 
wave source by a plurality sound ray vectors; 

storing propagation history data of a sound ray vector which data relates to one of said boundaries which is 
within a distance that said sound waves travel in a predetermined period of time and upon which said sound 
ray vector is incident to become an incident sound ray vector and at which said sound ray vector is reflected 
to become a reflected sound ray vector, wherein said propagation history data comprises said incident sound 
ray vector, said reflected sound ray vector, a total propagation distance between said wave source and said 
one of said boundaries, and coordinates of said intersection at which the sound ray vector intersects with said 
35 one of said boundaries; and 

determining said acoustic characteristics on the basis of both said stored propagation history data and a micro- 
area of said one of said boundaries occupied by said sound ray vector corresponding to said propagation 
history data. 
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9. A method as claimed in Claim 8, wherein said acoustic characteristics of said sound affecting said virtual window 
at predetermined time intervals are added to a time-series numerical array corresponding to said predetermined 
time intervals and stored, so that a transient response of said sound appearing at said observation point is deter- 
mined. 

10. A method as claimed in Claim 9, wherein a plurality of loud speakers are disposed in positions corresponding to 
said areas of virtual window to reproduce said stored transient response of said sound as to each of aid virtual 
window's areas, whereby sound field synthesis at said observation point is realized. 

11. Apparatus for reproducing acoustic characteristics of a three-dimensional virtual space sound at an observation 
point in a virtual space in which sound waves issued from at lease one sound source are propagated to reach said 
observation point, characterized by a first memory area for storing the number of a plurality of boundaries of a 
polygon entirely defining said virtual space, the coordinates of each of said boundaries, the coordinates of an 
intersection at which, said boundaries intersect with each other, and a plurality of sound ray vectors radiating from 
at least one wave source: 

a first processing means for determining the propagation history data comprising an incident sound ray vector, 
a reflected sound ray vector, the total propagation distance between said wave source and one of said bound- 
aries intersecting with said sound ray vector, and coordinates of an intersection at which said sound ray vector 
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intersects with said one of boundaries, said one being within a distance that said sound waves travel in a 
predetermined period of time, upon which one said sound ray vector is incident to become said incident sound 
ray vector, and at which one said sound ray vector is reflected to become said reflected sound ray vector: 
a second memory area for storing said propagation history data for said one of said boundaries; and 
a second processing means for determining said acoustic characteristics of said sound appearing at said 
observation point on the basis of both said stored propagation history data and a micro-area of said one of 
the boundaries occupied by said sound ray vector corresponding to said propagation history data. 

2. Apparatus as claimed in Claim 1 1 ; wherein said first memory area is arranged to store coordinates of at least one 
virtual window divided into a plurality of areas, and said virtual window from which said observation point is oppo- 
sitely-disposed in said virtual space. 

3. Apparatus as claimed in Claim 11 or 12, and further comprising a memory means; and 

said acoustic characteristics of said sound affecting said observation point at predetermined time intervals are 
added to a time-series numerical array corresponding to said predetermined time intervals and stored in said 
memory means; 

whereby a transient response of said sound appearing at said observation point is stored in said memory 
means. 

4. Apparatus for reproducing the three-dimensional virtual space sound, as claimed in Claim 13, wherein the appa- 
ratus further comprises a reproduction means in which a plurality of loud speakers are disposed in positions cor- 
responding to said areas of virtual window to reproduce by each of said corresponding loud speakers said transient 
response of said sound stored as to each of said virtual window's areas, whereby sound field synthesis at said 
observation point is realized. 

5. Apparatus as claimed in Claim 13, wherein a plurality of loud speakers for reproducing said virtual space sound 
in a real space are arranged in the apparatus; 

a plurality of sound directions in which said sound waves reach a listener are defined in the apparatus in 
accordance with positions of said loud speakers; and 

said transient response of said sound is determined for each of said plurality of sound directions. 

6. Apparatus for reproducing a three-dimensional virtual space sound, as claimed in Claim 15, wherein said transient 
response of said sound is determined for each of the combinations of said at least one sound source and said 
sound directions. 

7. Apparatus as claimed in Claim 16, and further comprising a sum-of-products calculator for reproducing said tran- 
sient response of said sound thus determined for each of said combinations, whereby a sound field for said ob- 
servation point is reproduced in said real space. 

8. Apparatus for reproducing a three-dimensional virtual space sound, wherein sound waves are issued from at least 
one sound source and propagated in a virtual space defined by a plurality of boundaries so that said sound is 
reproduced at an arbitrary observation point in said virtual space, characterized in that at least one virtual window 
is oppositely disposed from said observation point and which virtual window is divided into a plurality of areas; and 
wherein acoustic characteristics of said sound passing through each of said areas of virtual window are determined. 

9. Apparatus as claimed in Claim 1 8, further comprising a plurality of loud speakers disposed in positions correspond- 
ing to said areas of virtual window to reproduce said acoustic characteristics of said sound being determined for 
each of said areas of virtual window, whereby sound field synthesis at said observation point is realized. 
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processed exceed N ? 

Was the transient response calculated 
at every observation point ? 
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FIG.4 



Start 



Calculate the First integration term of equation 3 



41 



42 



43 



Calculate the second integration term of equation 3 



Calculate, at the point Q of the boundary B, an 
area A S of the boundary surface represented 
by the sound ray vector 



44 



Determine the product of the calculation result 
of step 41 and the A S , and process the product 
and the initial value of the sound sourceby the use 
of a convolution transform 



Determine the product of the calculation result 
of step 42 and the A S , and process the product 
and the derivative value of the sound source's 
initial value by the use of a convolution transform 
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46 



Calculate the arrival time Dt of the wave motion 
to reach the observation point P, based on 
distance DD which is the sum of the total 
propagation distance d up to the intersection Q 
and the distanrp F?D 



47 



Determine the array position J corresponding to 
the arrival time Dt in a numerical array storing 
the transient response that is being calculated 



48 



Sum up and store the time series data from 
steps 44, 45 in the array position J in which the 
transient response is stored; and, when some value 
Ls already stored in the array position J and the 
like, sum up the time series data and such value 



End 
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FIG.5A 




Potential issued from 
the micro-area element to 
afFect the observation point 



FIG.5B 



An area size of the micro-area element 
in the wave front, defined by a single 
sound ray vector 



Sound ray vectors 



Micro-area element 



Boundary 




Reflected sound ray vector 



o 

Observation point 



Potential issued from 
the micro-area element to 
affect the observation point 
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FIG. 6 



Sound ray vector 1 



Sound ray vector 2 



Sound ray vector 3 



Sound ray vector N 



Transient response 
to be determined 



Potential calculated at the observation point 
each time the sound ray vector is reflected at 
the boundary 



v 



^ — di t- 



rff ^V 0 . 




Simultaneous potentials 
are summed up 



th 



Direct sound R e nected sound 
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FIG.7A 



Wave sourc 




Potential from the boundary 

Listener 



FIG.7B 

Potential from the boundary 




Listener 



FIG.7C 



' Acoustic generator (loud speaker) 

o'o, 



0 Q<?/o 0 u ° 



0 



o P 0 



0 



0 



Sound waves 
.synthetically produced 



o o 88 oggr i 



Listener 



A large number of acoustic generators so arrauged 
as to surround the listener in any direction 
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FIG. 8 



Sound source such as 




M pieces of impulse responses having been 
calculated are registered in the sum of products 
calculator 
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FIG.9A 



Virtual space formed 
in the computer 

Sound waves from the virtual space 




Virtual window 



A room in real space 



Example of a virtual window sandwiched between 
virtual space and real space 



FIG.9B 



Virtual window constructed 
of a large number of the loud 
speakers 



Sound propagated from 
the virtual space 




Virtual space formed 
in the computer 



Real 



space 
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FIG. 10 




Example in which the loud speakers forming the 

virtual windows are hidden by a screen permeable 

to sound. The screen may serve as a picture-display screen 
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FIG.11A 




Virtual window constructed of 96 pieces of the 
loud speakers serving as the acoustic generators 



FIG.11B 




of products calculator 
Sura of products calculator with 24 channels 
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FIG. 13 

Acoustic generator 
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FIG. 14 

Wove 
source^ 

Sound ray vector 



Sound ray vector 




Sound ray vector 



^-Potential from the 
rnicroarea on the 
boundary surface 



' Sound ray vector 



— Observation 
point 



^Boundary of the 
virtual sound field 



Concept view of the sound field, 
taken from above 



FIG. 15 



Micro-area on the boundary 



Potential from the micro-area 
on the boundary surface 



Acoustic generator 1 



position of the acoustic 
generator (loud speaker) 
in the real sound field 



Acoustic generator 4 % 

1 4 Y 




'1x7;- Acoustic generator 2 



position of the acoustic 
generator (loud speaker) 
in the real sound field 



Observation point 
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FIG. 16 




Listener 
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FIG. 17 



Micro-area on the boundary 




Micro-area on the boundary 




Determine a horizontal component of 
the potential " <t>" from the boundary 
surface by calculating a function of 
0cos0 and integrating the result 
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FIG. 19 



Micro-area on the boundary 




Potential from the micro-area 
of the boundary surface 



Loud speaker 1 
Observation point 



Loud speaker 3 
Observation point 4 




Loud speaker 2 
Observation point 2 



Central part of the listener's ar 

^ Loud speaker 4 
"Observation point 3 
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FIG. 21 



Design the virtual space which a user wants to 
enjoy through synthesis of the sound field: imaged 
space, actually existing space such as concert 
halls and the like 



Numerically represent the virtual space thus 
designed 



Determine hardware to be used in sound field 
synthesis 



By the use of an approximate calculation of 
KirchhofTs integral, calculate the transient 
response of the propagated sound waves at the 
observation point in the virtual space 



By the use of both, the sound field synthesis 
hardware and the transient response thus 
calculated, synthetically produce a synthetic 
sound field of the virtual space in real space 
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55 
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FIG. 22 



Reflected sound 




Sound source 



Reflected sound 



(Elevation view) 



A numerically represented virtual space such 
concert halls and the like 



as 



^^ReQected sound 




Sound source 



(Plan view) 



A numerically represented virtual space such as 
concert halls and the like 
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FIG. 23 
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FIG.25A 
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Response obtained using a conventional orthodox 
calculus (virtuaed image method) 



FIG.25B 

a ODD 
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time(sec) 



Impulse response obtained using the calculating 
procedure of the present invention 
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FIG.26 
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